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The nonclassicality of single photon-added thermal states in the thermal channel is investigated 
by exploring the volume of the negative part of the Wigner function. The Wigner functions become 
positive when the decay time exceeds a threshold value yt c , which only depends on the effective tem- 
perature or mean thermal photon number of the thermal channel, but not depends on the effective 
temperature of the initial thermal state. This phenomenon is similar with the case of single quantum 
excitation of classical coherent states in thermal channel. Furthermore, we firstly demonstrate 7i c 
is the same for arbitrary pure or mixed nonclassical optical fields with zero population in vacuum 
state. 
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The preparations of nonclassical non-gaussian optical 
fields have attracted much attention, which have many 
applications in quantum information processing d| . Usu- 
ally, the nonclassicality manifests itself in specific prop- 
erties of quantum statistics 0, d, 0, 11 , in which the par- 
tial negative Wigner distribution [5| is indeed a good 
indication of the highly nonclassical character of the 
optical fields. Reconstruction of the Wigner distribu- 
tion in experiments with quantum tomography 0, @, Q 
have demonstrated the appearance of the negative val- 
ues, which does not have any classical counterparts. A 
variety of nonclassical states has recently been character- 
ized by means of the negativity of their Wigner function 
0, [l^, HH, EH, EH > which is a sufficient and not necessary 
condition for nonclassicality ■ 

The photon-added thermal states (PATSs) were intro- 
duced by Agarwal and Tara [i~5l |. which do not exhibit 
squeezing, sub-Poissonian statistics, and any coherence. 
The single photon-added thermal state (SPATS) has been 
experimentally prepared by Zavatta et al. and its non- 
classical properties have been detected by homodyne to- 
mography technology [l6j]. For the SPATSs, their non- 
classical properties have been investigated by several au- 
thors [HI, EH E3, EH- Furthermore, Parigi et al. have 
experimentally investigated quantum commutation rules 
by addition and subtraction of single photons to or from 
a light field initially in the thermal state [l9( . Ordinarily, 
the interaction between the nonclassical optical fields and 
their surrounding thermal environment may deteriorate 
the degree of nonclassicality. Thus, to study the dynam- 
ical behaviors of the partial negativity of Wigner distri- 
bution and understand how long a nonclassical field pre- 
serves its partial negative Wigner distribution in thermal 
channel may be very desirable for experimentally quanti- 
fying the variation of nonclassicality. Here, the nonclas- 
sicality of photon-added thermal states in the thermal 
channel is investigated by exploring the partial negative 
Wigner distribution. The exact expression of the time 
evolution Wigner function is given out and the dynam- 
ical behavior of the volume of the negative part of the 
Wigner function is analytically derived. The threshold 



values jt of the decay time corresponding to the tran- 
sition of the Wigner distribution from partial negative 
to completely positive are derived. For SPATSs in ther- 
mal channel, it is shown that the threshold value of the 
decay time is independent of the mean thermal photon 
number of initial thermal state seed. Furthermore, we 
firstly demonstrate "ft c is the same for arbitrary pure or 
mixed nonclassical optical fields with zero population in 
vacuum state, and is given by 7i c = m (|±|^), where n is 
the mean thermal photon number of the thermal channel. 

Let us first briefly recall the definition of the sin- 
gle photon-added thermal states (SPATSs) The 
SPATSs are defined by 
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where \l) is the Fock state and n is the mean photon num- 
ber of the thermal state seed. When the SPATS evolves 
in the thermal channel, the evolution of the density ma- 
trix can be described by [2(| 
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where 7 represents dissipative coefficient and n denotes 
the mean thermal photon number of the thermal channel. 
a^ (a) is the creation operator (annihilation operator) of 
the optical mode. When n = 0, the Eq.(2) reduces to the 
master equation describing the photon-loss channel. 

For an optical field in the state p, its Wigner function, 
the Fourier transformation of characteristics function [2l| 
of the state p can be derived by [22j, |23| 
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where 6 e = J2T=o |2fc)(2fc| and 6 Q = J2kLo |2fc + l)(2fc + 
1| are the even and odd parity operators respectively. 
In the thermal channel described by the master Eq.(2), 
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the time evolution Wigner function satisfies the following 
Fokker-Planck equation [24] 
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where and p represent the real part and imaginary part 
of /3, respectively. The time evolution Wigner function 
can be derived as following: 
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where 
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is the Wigner function of the thermal state with mean 
photon number n. Substituting the initial Wigner func- 
tion of the SPATS Ft! 
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into the Eq.(5), it is easy to obtain the corresponding 
time evolution Wigner function as follows: 
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In Fig.l, the Wigner function of the SPATS with n — 
1 in the thermal channel with n — 0.5 is plotted for 
three different values of decay time. It is shown that the 
phase space Wigner distribution of the SPATS exhibits 
partial negativity around the origin, and the region of 
the negative part in phase space is a circle. The absolute 
value of negative minimum of the Wigner distribution 
decreases as jt increases, and the thermal noise causes 
the disappearance of the partial negativity of the Wigner 
function if the decay time exceeds a threshold value. The 
ringlike wings in the distribution gradually disappear and 
the distribution becomes more and more similar to the 
Gaussian typical of a thermal state. 

Recently, the volume Pnw of the negative part of 
Wigner distribution function has been suggested as a 
good choice for quantifying the nonclassicality [H, [26|, 
|2ZL IH HI - Pnw is defined by 
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FIG. 1: The Wigner functions of the SPATS with n = 1 in 
thermal channel with n — 0.5 are depicted for three different 
values of decay time jt. 



where f2 is the negative Wigner distribution region. In 
Ref. [HI , we have investigated Pnw 01 photon-added co- 
herent states in the photon-loss channel. It was shown 
that Pnw and entanglement potential defined in Ref. [30j 
exhibit the consistent behaviors in short decay time. 

Now, we bring our attention to the influence of ther- 
mal noise on the nonclassicality of the quantum excita- 
tion of classical non-coherent optical thermal fields by 
calculating Pnw ■ Substituting the time evolving Wigner 
function in Eq.(8) into Eq.(9), we could obtain 
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lation between the mean photon number n of the thermal 
channel and the threshold decay time r yt c (n) [29], 



FIG. 2: The volume of partial negativity of the evolving 
Wigner function of the SPATSs with different values of n in 
the thermal channels with two different values of mean ther- 
mal photon number (a) n = 0.5; (b) n — is plotted as the 
function of decay time 7 t. (Solid line) n = 0; (Dash line) 
n = 2; (Dot line) n = 1. 



for 7* < 7 i c = lnf±fg. In Fig.2, P ww of SPATSs with 
different values of the parameter n is plotted as the func- 
tion of "ft for two different values of the mean thermal 
photon number n of the thermal channel. It is shown that 
the thermal noise deteriorates the partial negativity, and 
Pnw monotonically decreases with the decay time. The 
more intense the thermal noise, the more rapidly Pnw 
decreases, which implies that the nonclassicality of the 
optical fields are very fragile against the thermal noise. 
The above results also indicate that Pnw becomes zero 
at a threshold decay time jt c which only depends on the 
value of n but not n. For the case of the SPATSs in 
thermal channel, the threshold decay time is 
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Previous studies have indicated that the volume of neg- 
ative part of Wigner function of single photon-added co- 
herent states (SPACSs) in photon-loss channel or thermal 
channel decreases with the decay time, and its thresh- 
old decay time concerning the complete vanish of partial 
negativity of Wig ner function does not depend on the 
seed intensity [28|, [29| . In Ref. [29J] , we have presented a 
relation between the mean photon number n of the ther- 
mal channel and the threshold decay time "/t c (n) beyond 
which their Wigner function become positive. For ar- 
bitrary nonclassical optical fields which have the partial 
negative Wigner distribution function, there exists a re- 
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where 7t c (0) is the threshold decay time in the photon 
loss channel. 

The similarity between the amplitude-independence of 
the threshold decay time 7 i c corresponding to the dis- 
appearance of partial negativity of the Wigner function 
of SPACSs and SPATSs in the thermal channel implies 
there may exist a universal relation about 7 i c for ar- 
bitrary single quantum excitation of classical gaussian 
states. Here, we demonstrate that all of the nonclassi- 
cal states whose density operators p satisfy (0|p|0) = 0, 
where |0) is the vacuum sate, completely loss the nega- 
tivity of their Wigner functions at the threshold decay 
time 7 i = In fjrfjj ■ The strict proof will be presented 
elsewhere. Here, we briefly outline the proof procedures. 
Firstly, considering the case of photon-loss channel, i.e. 
n = 0. In this situation, based on Eq.(5), for arbitrary 
initial nonclassical pure or mixed states whose vacuum 
state population is zero, we can derive 

7i— In 2 > o (13) 

for any values of q and p, and 

W(0,0,7t)| 7t=ln 2=c<0M0)=0 (14) 

where c is a constant. Eqs. (13-14) imply that the thresh- 
old decay time 7 £ c (0) in Eq.(12) is In 2 for these states. 
Substituting it into Eq.(12), we can complete this proof. 
In the derivation of Eq.(13), we have used the relation 
between the Wigner distribution function and the Q 
function. The Q function gives the probability distri- 
bution for finding the coherent state \a) in the state p 
since Q(q,p) = ~(a\p\a), where a = q + ip. The Q 
function is always non-negative. The Q function is a 
particular case of a class of non-negative quantum dis- 
tribution functions, the Husimi functions, obtained by 
smoothing the Wigner distribution function with a min- 
imum uncertainty squeezed Gaussian function 31]. The 
Q function can be obtained when the Wigner function is 
smoothed by a coherent state wave packet. From Eq.(5), 
we can find, for any quantum fields with initial state 
po in photon loss channel, their time evolution Wigner 
functions at decay time yt = In 2 can be rewritten as 
W(g,p,ln2) cx (V2a\p \V2a) = 7rQ (V2q, y/2p). There- 
fore, we have the Eq.(13) by referencing the characteris- 
tics of the Q function. 

In summary, we have investigated the nonclassicality 
of single photon excitation of thermal optical field in the 
thermal channel by exploring the partial negativity of 
the Wigner function. The total volume of the negative 
part defined by the absolute value of the integral of the 
Wigner function over the negative distribution region is 
analytically calculated. For the case of SPATSs in ther- 
mal channel, the exact threshold value of the decay time 
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beyond which the evolving Wigner function becomes pos- 
itive is given as jt c = ln( ^+^" ), which is the same as the 
one in the case of single quantum excitation of the clas- 
sical coherent field. 

For all of the nonclassical states whose density opera- 
tors p satisfy (0|p|0) = 0, where |0) is the vacuum sate, 
it is demonstrated that the threshold decay times are 
the same and given by 7 i c = ln(|±|n). These results 
clearly imply Wigner distributions of any photon-added 
thermal states are partial negative before the threshold 
decay time "ft c even if their initial thermal state seeds 
are macroscopic with arbitrary large n but finite. Obvi- 
ously, any photon-added classical gaussian states includ- 
ing photon-added thermal states, photon-added coherent 
states, and photon-added displaced thermal states belong 
to the class of states satisfying (0 1 p| 0) = 0. Therefore, 
the above conclusions can be generalized to: In thermal 



channel with mean thermal photon number n, all non- 
classical pure or mixed states p with zero population in 
vacuum state have partial negative Wigner distribution 
before the threshold decay time jt c = In jjf^ if Tr(pa'a) 
is finite. 

The above results can be regarded as benchmark to 
investigate the robustness of other indicators of nonclas- 
sicality such as squeezing, antibunch, and entanglement 
potential of nonclassical optical fields in thermal chan- 
nel if compared with the partial negative Wigner distri- 
bution. Recently, the physical realization of controlled 
phase gate based on the single-photon-added coherent 
states has also been proposed [29|]. Our studies in this 
report may find some applications in these quantum in- 
formation processes in which photon-added states are in- 
volved. 
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